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ABSTRACT 

Research at the Naval Postgraduate School has led to the develop- 
ment of a system for producing GPSS simulation programs for simple 
queuing problems through English language dialogue with an IBM 
360/67 computer. This thesis describes work done to give the system 
the capability to actually perform the simulation and report the results 
through English language dialogue. A complete sample terminal 
session is included. 
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I. INTRODUCTION 



Since the advent of the computer, men have been striving to find 
easier methods of communicating their problems to the machine. 
Primitive communication was established by using the language of the 
machine and conversing at the computer’s elementary level. In an 
effort to narrow the communication gap between man and machine, 
translators have been developed. These translators - assemblers and 
compilers - facilitate communication with the computer at a level 
considerably removed from the machined language. However, even 
at this level, man is required to learn a new language in order to 
converse with his powerful assistant. Although significant advances 
have been made toward generalized, multi-purpose, higher level 
languages, considerable effort must be expended to learn and utilize 
these languages in a problem solving situation. A desirable alterna- 
tive would be to have the ability to specify the problem directly to the 
machine in a natural language such as English, and have the machine 
solve the problem and report the results as requested by the user. 

Several research projects have investigated various aspects of 
natural language interaction between man and computer [l, 2]. Develop- 
ments in the fields of linguistics and artificial intelligence have served 
to provide basic conceptual structures for natural language processing. 
One such development is the theory of stratificational linguistics 
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proposed by S. M. Lamb [3]. In this theory, language is considered 
to be a multi-level system of relationships. 

A research project is currently being conducted at the Naval Post- 
graduate School to investigate using natural language man-machine 
interaction in solving queuing problems by simulation [4]. The system 
being developed is called NLPQ and is a specific application of a more 
general system called NLP. This work is based on the concepts of 
s tratificational linguistics and utilizes an entity-attribute-value data 
structure. Background information on the development of both systems 
is included in this chapter. A further description of each system is 
included in Chapter II. 

A. BACKGROUND 

The initial work done on the project was the development of the 
general system NLP or Natural Language Processor. The constituents 
of this system are a "rule language" and a set of FORTRAN routines 
which compile and execute statements in the rule language. System 
monitor functions are also performed by the main routine. The system 
is implemented on the IBM 360/67 and is executed under control of the 
CP/CMS time sharing system. 

With the basic system established, further research began to 
produce the rule modules necessary to handle a queuing system 
application (NLPQ). The form of an internal problem description (IPD) 
for a queuing problem was decided upon and a set of encoding rules 
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were written to convert the information contained in an IPD to a GPSS 
program [5]. Additional encoding rules were added to produce an 
English text description of the information contained in an IPD [6]. 

The generality of the English encoding rules also permits their use 
in other areas involving natural language responses from the computer. 
A set of English decoding rules was then developed to allow the user 
to describe his queuing problem to the computer in English. These 
rules perform the function of processing the input English text to 
produce the IPD. In addition, they are utilized in handling English 
language requests from the user. 

Further research developed modules which would massage and 
inspect the IPD for missing or erroneous information before producing 
a GPSS program [7, 8]. In cases where missing or erroneous informa- 
tion is detected, a request is made to the user to supply or correct 
the required information. A practical by-product of this research is 
the capability to enter an English problem description in a question- 
answer mode. 

More recently, a FORTRAN subroutine designed to perform a 
GPSS-like simulation and a set of encoding rules to initialize the data 
structure required by this routine were developed [9, 10]. Information 
contained in the IPD can be manipulated by these rules to produce 
either a GPSS program or a representation of the queuing problem in 
the data structure utilized by the simulation routine, or both. 
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B. THESIS OBJECTIVE 



The objective of the research for this thesis was to integrate the 
simulation routine and associated rules into the existing NLPQ system 
to produce an initial version of an interactive simulation system. 

This required making modifications and extensions to both the 
simulation routine and several of the existing rule modules. 

C. ORGANIZATION OF THE THESIS 

Chapter II of this thesis presents a more detailed discussion of 
pertinent portions of NLP and NLPQ. Chapter III contains a sample 
session to illustrate the capabilities of NLPQ in an interactive problem 
solving situation. The considerations involved in the implementation 
of the interactive simulation capability are discussed in detail in 
Chapter IV. Finally, Chapter V presents conclusions and recommenda- 
tions for further research. 
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II. DESCRIPTION OF NLP AND NLPQ 



Since the interactive simulation capabilities are integrated with, 
and rely on, the other components of NLPQ, a description of those 
modules will be presented in this chapter. First, however, the basic 
concepts related to an overview of the general system NLP, the data 
structure utilized, and the rule language will be discussed. 

A. BASIC CONCEPTS 

This section is intended to provide a general outline of the basic 
concepts inherent in NLP and NLPQ. A detailed discussion of this 
material may be found in Ref. 4. 

NLP is composed of a set of FORTRAN routines and a rule lan- 
guage. The main program serves as a monitor and performs certain 
input/output operations. Subroutines compile statements in the rule 
language and interpret operations given by the rule statements. An 
entity-attribute- value data structure is utilized to hold information. 
The generality and usefulness of this type of data structure has been 
widely recognized in the fields of simulation programming systems 
and artificial intelligence. 

The entity-attribute-value data structure is well suited for holding 
several types of information. For example, in the current queuing 
problem application, information about the various words and concepts 
related to queuing problems must be maintained. Relations between 
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words and concepts can be considered to be held in long term memory, 
since information of this type is necessary to carry on a dialogue 
about any problem. 

Other information about a specific problem being described must 
be retained for the duration of the problem solving session. Information 
of this type is obtained through discourse. The problem description 
input by the user is first processed by the ’’decoding 11 rules. These 
rules serve to convert the information contained in the description to 
an equivalent internal representation. This type of storage can be 
considered as short term memory since the system need only retain 
it for the duration of the specific problem solving session. 

Another type of information storage can be considered to be 
temporary or scratchpad memory. For instance, information about 
parts of a sentence can be discarded after the sentence has been 
completely processed. An important aspect of NLP is that all of these 
types of information are maintained in exactly the same form, i. e. , 
entity-attribute- value. Thus, all types of information can be manip- 
ulated in the same fashion. 

Rules in the rule language of NLP consist of two parts separated 
by an arrow. The left part generally specifies conditions which must 
be satisfied before the rule can be applied. The right part of the rule 
specifies the actions to be taken when the rule is executed. Various 
basic elements, known as records, establish the state of the system. 
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These records carry the types of information mentioned in the pre- 
ceding paragraphs. 

Those records which are used in a scratchpad fashion to create 
or modify other records are called segment records. The state con- 
ditions contained in this type of record are the most frequently tested 
by the rules. The rules for "decoding 11 specify the manner in which 
records are to be created from input character strings. The "encod- 
ing" rules describe the inverse conversion from records to character 
strings. Record-to-record transformations can be accomplished by 
either type of rule. Thus, the basic system deals with the conversion 
of information. Information in the form of a natural language character 
string may be transformed to an internal format, manipulated, and 
possibly returned to a character string. The modular sets of rules 
for performing these functions for NLPQ will be considered below. 

Since the internal problem description (IPD) is the structure to be 
built by the decoding rules, it will be discussed first. 

B. THE INTERNAL PROBLEM DESCRIPTION 

Utilizing the entity-attribute -value data structure previously 
discussed, the logical structure of the IPD is a set of records which 
contain information about the current problem being processed by 
NLPQ. These records represent entities such as physical objects or 
actions occurring in the problem. These entities have attributes which 
in turn have values associated with them. 
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A typical queuing problem sequence involves mobile entities 
engaging in actions (abstract entities) at various stationary entities. 

In most cases, each of these entities has attributes of varying com- 
plexity with specified values. Each of these records in the IPD is a 
member of one of seven lists. Actions are members of the action 
list ( 'ACTNLIST 1 ), mobile and stationary entities are members of the 
‘MOBLIST 1 and 'STALIST', respectively. The other lists are: the 
distribution list ( 'DSTRLIST' ), the successor descriptor list 
( 'SCSRLIST'), the miscellaneous list ( 'MISCLIST '), and the unit list 
( *UNITLIS T 1 ). Each of these lists is a "named record". Named 
records provide the long term memory capability previously described; 
that is, they contain word and concept information pertinent to the 
application. These particular list records, however, are used to hold 
information about a specific problem. As entities are encountered 
during decoding, records are created and linked into the appropriate 
lists. 

The concept structure created by the named records provides the 
framework of words and their semantic content necessary for discourse. 
As such they represent the system's vocabulary and knowledge of the 
relationships between words and concepts. Using this general knowl- 
edge, a "mental image" (the IPD) of the specific problem entered by 
the user can be obtained. The IPD can be considered to be a specific 
problem instance in the domain of the queuing conceptual structure. 
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C. DECODING THE PROBLEM DESCRIPTION 

In NLPQ the decoding rules specify how input text is to be con- 
verted into equivalent information in the form of records. This 
conversion is performed within the framework of the Stratificational 
Grammar theory. Within this theory several levels (strata) of language 
structure exist. Three levels have been utilized in the NLPQ applica- 
tion; the morphological, lexological, and the semological levels. The 
n morphology !! is concerned with the manner in which characters are 
put together to form words. As such, it is highly dependent on the 
particular language being used. The n lexology n deals with the way in 
which words are put together to form phrases, clauses, and sentences. 
Different grammatical orderings required by different languages 
result in language dependency at this level also. The semological 
level, however, is concerned with relationships and meanings. Thus 
elements at this level are relatively language independent. The decod- 
ing process then involves applying morphological and lexological rules 
first in processing the input text. Semological rules are then utilized 
to develop the structure representing the meaning of the text, the IPD. 

The general format of a decoding rule is: 

SEGMENT TYPE (COND 1, 2, . . . ) SEGMENT TYPE (COND 1, 2, ...)... 

^ SEGMENT TYPE (ACTN 1, 2, ... ) 

Thus a decoding rule specifies what to do when a particular series of 
segment types (satisfying the given conditions) is found while processing 
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the input text. Rule application results in the creation of the segment 
type on the right and performance of those actions indicated. The 
conditions specified on the left side of the rule are known as "con- 
dition specifications". The actions performed in creation of the new 
segment on the right side are known as "creation specifications". 

Both types of specification elements have access to and can manipulate 
any information in the system. 

i 

An illustration of some of these features can be seen in the follow- 
ing rule from the decoding morphology: 

VERBS(ING) I N G VERBP(SUP( VERBS), PRESPART) 

The left part of the rule is made up of four segment types, a verb 
stem (VERBS) segment and three "character" segments. The con- 
dition specification for the verb stem requires that the ING indicator 
in that segment be "on". Indicators are binary-valued and indicate 
the presence or absence of certain attributes. The right part of the 
rule defines the conversion to be made when a series of segment types 
satisfying the left part is encountered. In this case, a new verb part 
(VERBP) segment is created which is in the same superset as the 
verb stem and has a present participle indicator on. 

Using this basic scheme, information is extracted from the input 
text and used to build the IPD. Once the semantic content of the user's 
dialogue has been transformed to the internal format, it can be manip- 
ulated in several ways. 
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D. ENCODING FROM THE IPD 



Since encoding is basically the inverse process of decoding, 
encoding rules are essentially the inverse of decoding rules. Informa- 
tion is manipulated from the semological level, through the lexological 
and morphological strata to a natural language text output. A general- 
ized format for an encoding rule is: 

SEGMENT TYPE (COND 1, 2, ... ) * ^ 

SEGMENT TYPE (ACTN 1, 2, ... ) SEGMENT TYPE (ACTN 1,2,...) . . . 
In this case, the rule specifies the sequence of segment types (with 
corresponding actions) to be created when a segment type satisfying 
the appropriate condition specifications is encountered. 

Both encoding and decoding rules have the ability to perform 
record-to- record information conversion, as previously stated. One 
modular set of encoding rules known as the MASSAGER [7] is utilized 
in this way to set default values in the IPD. Certain assumptions about 
the problem may be made by the user during the discourse. For 
example, if the number of units of storage capacity required by a 
mobile entity has not been mentioned, it is probable that an assumption 
of one unit has been made by the user. The purpose of the MASSAGER 
is to inspect the IPD after decoding is completed and set default values 
in those instances where non- controversial assumptions can be made. 
Other functions include the consolidation of redundant information and 
the deletion of certain attributes required only for decoding purposes. 
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Another encoding rule module known as the INTERROGATOR [8] 
serves to inspect the IPD for missing or erroneous information. 
Copies of the action records maintained in the IPD are accessed 
through the 'ACTNLIST 1 and tested to ensure that they have all of the 
attributes required for processing by the GPSS /X- VECTOR rule 
module. When incorrect or missing information is detected, the 
INTERROGATOR sets up the segment form of the question to be asked 
and invokes the ENGLISH encoding module to actually produce the 
question. Information provided by the user is then decoded and the 
inspection of the IPD continues. The INTERROGATOR may also be 
utilized to enter the problem in a question-answer fashion. When the 
necessary information pertaining to the problem has been established, 
a message is encoded indicating completion of the problem statement. 

At any point during the discourse the user may request a state- 
ment of the problem as the system "sees" it. The English encoding 
module [6] provides the conversion necessary to produce an English 
description of the problem from the information contained in the IPD. 
The generality of this module also permits the handling of statements 
to the user generated by other modules such as the INTERROGATOR. 

The GPSS/X- VECTOR encoding module [ 1 0] is utilized to create 
a GPSS program and initialize the data structure used by the simula- 
tion routine. This data structure, called the X-vector, contains the 
information necessary to perform the simulation. In addition, it is 
used throughout the simulation to maintain the required statistics in 
much the same manner as the internal tables of GPSS [ii] 



17 



Once the problem has been completely specified, the user may 
request that a GPSS program be written for the current problem. The 
semological rules of this module examine the IPD and produce segments 
which roughly correspond to the statements of a GPSS program. Rules 
in the lexology further process these segments and their constituents 
into other segments in the appropriate order for a GPSS program. 

The morphological rules then produce the corresponding GPSS output. 

Similar rules in the X-vector lexology and morphology place 
pertinent information into the X-vector. Figure 1 (taken from ref. 10) 
shows an initial X-vector produced by these rules. With the informa- 
tion in the X-vector the simulation can be performed. 

E. THE SIMULATION ROUTINE 

This FORTRAN routine, originated by Williams [9], performs a 
GPSS-like simulation based on the information contained in the X-vector 
The initial portion of the X-vector contains "parameters" for the sim- 
ulation routine. These parameters are assigned fixed locations in the 
vector. They consist of variables such as the random number seeds 
to be used, pointers to the various directories, entity allocation counts, 
clock time, etc. The latter portion of the vector contains allocated 
space for the various GPSS entities (i. e. , STORAGES, QUEUES, 
TABLES, FUNCTIONS, VARIABLES, SAVEVALUES, and BLOCKS). 

The directories associated with these entities are also allocated 
space in this section of the vector. The location of these allocated 
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1 


MODE 
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2 


TERMINATION 

COUNT 


1 


3 


SEEDl 


277 


4 


SEED2 


423 


5 


SEED3 


815 


6 


SEED4 


121 


7 


SEED5 


655 


8 


SEED6 


531 


9 


SEED7 


999 


10 


SEEDS 


813 


11 


CLOCK TIME 


0 


12 


NUMBER OF 
STORAGES 


2 


13 


STORAGE 
DIR. PTR. 


336 


14 


NUMBER OF 
QUEUES 


2 


15 


QUEUE 
DIR. PTR. 


338 


16 


NUMBER OF 
FUNCTIONS 


4 


17 


FUNCTION 
DIR. PTR. 


343 


18 


NUMBER OF 
BLOCKS 


14 


19 


BLOCK 
DIR. PTR. 


348 


20 


NUMBER OF 
TABLES 


3 


21 


TABLE 
DIR. PTR. 


340 


22 


NUMBER OF 
SAVE VALUES 


0 


23 


SAVEVALUE 
DIR. PTR. 


0 


24 


NUMBER OF 
PARAMETERS 


3 


25 


TRANSACTION 
POINTER 362 


26 


FIRST TRANS . 
CEC. PTR. 0 


27 


LAST TRANS 
CEC. PTR. 


0 


28 


FIRST TRANS. 
FEC. PTR. 0 



29 

30 

31 

32 

33 



43 



51 



61 



69 



80 



91 



169 



262 



LAST TRANS. 
FEC. PTR. 0 


282 


ERROR CODE 0 


288 


NUMBER OF 
VARIABLES 1 


VARIABLE 
DIR. PTR. 347 




STORAGE 

ALLOCATION 

(STAT1) 


337 


QUEUE 

ALLOCATION 

(STAT1) 


339 


STORAGE 
ALLOCATION 
(PUMP 2) 


341 


QUEUE 

ALLOCATION 
(PUMP 2) 


344 


TABLE 2 
ALLOCATION 


348 


TABLE 3 
ALLOCATION 


349 


EXPON 

FUNCTION- 

ALLOCATION 


363 




NORMAL 

FUNCTION 

ALLOCATIONS 




ADDITIONAL 

FUNCTION 

ALLOCATIONS 





VARIABLE 

ALLOCATION 



BLOCK 

ALLOCATIONS 



STORAGE 

DIRECTORY 



BLOCK 

DIRECTORY 



32 



50 



QUEUE 

DIRECTORY 



42 



60 



TABLE 

DIRECTORY 



0 

68 

79 



FUNCTION' 

DIRECTORY 



VARIABLE 

DIRECTORY 



TRANSACTION 

ALLOCATIONS 



Figure 1: Internal Structure of the X-Vector 
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areas is flexible and is determined by the requirements of the 
specific problem. 

The procedure used to execute the simulation is essentially the 
same as that of GPSS. Raw results of the simulation, such as the 
cumulative time integrals for storages and queues, are stored in the 
X-vector elements associated with these entities. Several simulation 
modes comparable to the GPSS control cards SIMULATE/START, 
RESET, CLEAR, and START perform similar functions in the 
simulation routine. 
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III. A SAMPLE SIMULATION PROBLEM 



This chapter is intended to demonstrate the capabilities of this 
initial version of the interactive simulation system being developed. 

It illustrates how a user may input a queuing problem in English, 
have the system question him for needed information, obtain a restate- 
ment of the problem from the computer, and have the system perform 
the simulation and output the requested results. 

The sample given here is taken from an actual terminal session. 

In line with the purpose of demonstrating the system's capabilities, 
a wide variety of statements and questions are given. This results 
in an unusual amount of redundancy. 

Throughout this chapter, all inputs by the user are shown in 
lower case, and all computer responses are shown in UPPER CASE. 
This is the way a sample session under CP/CMS would actually appear. 
Comments are interspersed where needed and are enclosed in paren- 
theses. The sample session begins with "decode : 11 below. 



decode: 

(This is a command to NLP necessary to invoke the 
decoding rules for processing the input text. ) 

vehicles arrive at a station, 

the station has just one pump. 
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a vehicle will leave the station immediately after 
arriving if the length of the line at the pump 
is not less than two. 

otherwise, it is serviced there; then it leaves. 

service times are exponential, with a mean of 5 

minutes for cars and 9 minutes for trucks. 

three quarters of the vehicles are cars and one 

fourth of them are trucks. 

is the statement of the problem complete? 

(This question results in inspection of the IPD for 
missing or erroneous information. ) 

HOW OFTEN DO THE VEHICLES ARRIVE AT THE STATION? 

(The question-answer mode is initiated to query the 
user for necessary information. ) 

arrivals are normally distributed with a mean of 8 

minutes and a standard deviation of 2 minutes. 

HOW LONG SHALL THE SIMULATION BE RUN? 

the basic time unit to be used in the model is 

30 seconds. 

(Even though this response does not answer the question, 
this information is still accepted into the IPD. The 
system will query the user again for the answer to its 
question. ) 

HOW LONG SHALL THE SIMULATION BE RUN? 
eight hours. 

THE PROBLEM STATEMENT IS COMPLETE. 
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describe the problem in english. 



(This is requesting a restatement of the problem as it 
now appears in the IPD. ) 

THE VEHICLES ARRIVE AT THE STATION. THE TIME 
BETWEEN ARRIVALS OF THE VEHICLES AT THE STATION IS 
NORMALLY DISTRIBUTED, WITH A MEAN OF 8 MINUTES AND A 
STANDARD DEVIATION OF 2 MINUTES. 75 PERCENT OF THE 
VEHICLES ARE CARS, AND THE REST ARE TRUCKS. AFTER 
ARRIVING AT THE STATION, IF THE LENGTH OF THE LINE AT 
THE PUMP IN THE STATION IS LESS THAN 2, THE VEHICLES 
WILL BE SERVICED AT THE PUMP IN THE STATION. OTHERWISE, 
THE VEHICLES WILL LEAVE THE STATION. THE TIME FOR THE 
VEHICLES TO BE SERVICED AT THE PUMP IN THE STATION IS 
EXPONENTIALLY DISTRIBUTED, WITH A MEAN OF 5 MINUTES 
FOR THE CARS, AND 9 MINUTES FOR THE TRUCKS. AFTER 
BEING SERVICED AT THE PUMP IN THE STATION, THE VEHICLES 
LEAVE THE STATION. 

THE SIMULATION IS TO BE RUN FOR 8 HOURS, USING A 
BASIC TIME UNIT OF 30 SECONDS. 



write a gpss program for this problem. 

(Production of the GPSS program also results in X-vector 
initialization in preparation for running the simulation. 
The program produced is shown on the following page. ) 
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STAT1 

PUMP2 

CAR2 

2 

TRUC3 

3 

1 

0 . 0 , 0 . 

0.400, 

0.750, 

0.900, 

0.900, 

0.995, 



SIMULATE 

RMULT 277,423,715,121,055,531,999,813 

EQU 1/F,0 

EQU 2, F , Q 

EQU 2, T 

TABLE Ml, 1,1, 2 

EQU 3, T 

TABLE Ml, 1, 1 , 2 

FUNCTION RN1, C24 



0/0.100,0.104/0.200 

0.509/0.500,0.090/0 

1.390/0.800,1.000/0 

2.300/0.920,2.520/0 

3.200/0.970,3.500/0 

5.300/0.998,0.200/0 



,0.222/0.300 

.000,0.915/0 

.840,1.830/0 

.940,2.810/0 

.980,5.900/0 

.999,7.000/1 



0.355/ 

700,1.200/ 

880,2.120/ 

950,2.990/ 

990,4.000/ 

000 , 8 . 000 / 



2 FUNCTION RN2 , C29 

0.0,-3.000/0.012,-2.250/0.027,-1.930/0.043,-1.720/ 

0.002,-1.540/0.084,-1.380/0.104,-1.200/0.131,-1.120/ 

0.159,-1.000/0.187,-0.890/0.230,-0.740/0.207,-0.020/ 

0.334,-0.430/0.432,-0.170/0.500,0.0/0.508,0.170/ 

0.000,0.430/0.732,0.020/0.770,0.740/0.813,0.890/ 

0. 84 1,1.000/0. 809, 1.12 0/0. 890, 1.200/0. 910, 1.380/ 

0.938,1.540/0.957,1.720/0.973,1.930/0.988,2.250/ 

1.000,3.000/ 

3 FUNCTION PI, D2 

CAR2, 10/TRUC3, 18/ 

4 FUNCTION RM3 , D2 

0. 750, CAR2/1 . 000, TRUC3/ 

1 FVAR I ABLE 1G+4*FN2 



* THE VEHICLES ARRIVE AT THE STATION. 
GENERATE VI 

ASSIGN 1 , FN4 

TEST L Q$ PUMP2 , 2 , ACT 2 

TRANSFER ,ACT3 

* 

* THE VEHICLES LEAVE THE STATION. 

ACT2 TABULATE PI 
TERMINATE 



THE VEHICLES ARE SERVICED AT THE PUMP. 
ACT3 QUEUE PUKP2 

SEIZE PUMP 2 - 

DEPART PUMP 2 

ADVANCE FM3, FM1 

RELEASE PUMP 2 

TRANSFER ,ACT2 



* TIMING LOOP 

GENERATE 900 
TERMINATE 1 
START 1 

END 
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perform the simulation 



SIMULATION TIME IS 960(RELATIVE), 960(ABSOLUTE) . 



(This message signals completion of the simulation. The 
user may now ask for statistical printouts or for specific 
Information concerning the outcome. ) 

print the gpss statistics. 



STORAGE 


CAPACITY 


AVERAGE 

CONTENTS 


AVrRAGF 
UTI LIZATI0N 


FNTRI FS 


avtra nr 
T IMF/TRAN 


CURRENT 

CONTENTS 


MAXIMUM 

CONTENTS 


1 


1 


0.0 


0.0 


0 


0.0 


0 


0 


2 


1 


0.659 


0.659 


59 


10. 729 


0 


1 



QUEUE 


MAXIMUM 


A '/FRA OF 


TOTAL 


ZERO 


prRCENT 


AVFRACE 


SA'/FRAGF 


CURRENT 




CONTENTS 


CONTENTS 


E.UTRI fs 


ENTRIES 


ZEROS 


TIMF/TRANS 


TlMr/TPANS 


CONTENTS 


1 


0 


0.0 


0 


0 


0.0 


0. 0 


0.0 


o 


2 


2 


0.278 


59 


38 


64.407 


4.525 


12.734 


0 


$AVFRAGF 


timf/trans 


= AVFRACE 


T 1 NT /TRANS 


FXCl.Um fJO 


ZERO ENTRIES 









TABLE 2 



entries in tarlf 
43 



MFAN ARGUMENT 
13.628 



STANDARD nr vi ATI ON 
13.175 



SUM OF ARGUMENTS 
586.000 



UPPER 


OBSERVED 


PFR CFNT 


CUMULATI Vr 


CUMULATI VF 


LIMIT 


frequency 


OF TOTAL 


PFRCFNTACF 


REMAINDER 


1 


3 


6.98 


6.98 


93.02 


OVERFLOW 


40 


93.02 


100.00 


0. 00 


AVERAGE VALUF OF OVERFLOW 


14.650 



MULTI PLE 
OF MTAN 
0. 073 



DFVIATI0N 
FROM MEAN 
-0.958 



TABLE 3 












ENTRIES IN 


TAB L c 


MEAN 


ARGUMENT 


STANDARD DEVIATION 


18 






17.444 




12.641 


UPPER 


OBSERVED 


PFR 


: CENT 


CUMULATIVE 


CUMULATI Vr 


LIMIT 


FREQUENCY 


OF 


TOTAL 


PERCENTAGE 


Rr M A | f 'PFR 


1 


2 




11. 11 


11.11 


8 8.89 


OVERFLOW 


16 




88.89 


100.00 


o! 00 


AVFRAGF VALUF OF OVERFLOW 


19.625 







SUM OF ARGUMENTS 
314.000 



MULTI PLF 
OF MEAN 
0.057 



DFVI ATION 
FROM MFAN 
-1.301 



CURRENT FVFNTS CHAIN 



TRANS 


rot 


BLOCK 


NBA 


MARK-TIMF 


PI 


P2 


P3 P4 


SI 


Tl 


FUTURF 

TRAMS 


EVFNTS 

RPT 


CHAIN 

RLOCK 


NBA 


MARK-TIMF 


PI 


P2 


P3 P4 


SI 


Tl 


3 


964 


1 


2 


964 


0 


0 


0 23058912 


T 


T 


2 


1920 


13 


14 


1920 


0 


0 


0 372 


T 


T 
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print the storage statistics. 



STORAGE CAPACITY AVFRAOF AVrRAOF 

CONTENTS UTILIZATION 


FI.'TRI fs 


Avr ra r r 

T IMF/TRAN 


CURRENT 

CONTENTS 


MAXIMUM 

CONTENTS 


1 1 
2 1 


0.0 0.0 

0.659 0.659 


0 

59 


0.0 

10.729 


0 

0 


0 

1 


print the pump statistics. 










STORAGE CAPACITY AVFRAOF AVFRAOF 

CONTENTS IJTI LIZATION 


FMTRI FS 


AVFRAOF 

TIMF/TRAN 


CIJRRFNT 

CONTENTS 


MAXIMUM 

CONTENTS 


2 1 


0.659 0.659 


59 


10.729 


0 


1 


i 

OUFIIE MAXIMUM 

CONTENTS 


avfraof total zero 

CONTENTS FMTPIFS. ENTRIES 


PFRCTMT 

ZEROS 


A V C R A o r 

TIME/TRANS 


$ AVERAOF 
TIMF/TRANS 


CURRENT 

CONTENTS 


2 2 


0.278 59 38 


64.407 


4. 525 


12.714 


0 


JAVERAOF T IMF/TRANS 


- AVFRAOF TIMF/TRANS FXCLUPfNC 


ZFRO PMTRIES 








print the pump queue statistics. 










OUFUE MAXIMUM 

CONTENTS 


AVFRApc TOTAL ZFRO 

CONTENTS FNTP|r S FNTRIES 


pcpCFMT 
Z FRO S 


A V c R A G F 

TIME/TRAMS 


SAVERAOF 

TIMF/TRANS 


CURRENT 

CONTENTS 


2 2 


0.278 59 38 


64.407 


4. 525 


12.714 


0 



JAVERAOp TIMF/TRANS « AVFRApc T I (IT /TRANS FXC L'lR I f.'f! Z.TRO ^mtrIES 



print the truck statistics. 



TARLE 3 



ENTRIES III TARLF 
18 



MEAN ARCUMTUT STANOARR DEVI ATI ON 

17 * r, ' ,l ‘ 12.61.1 



UPPER OBSERVER 

LIMIT FREQUENCY 

1 2 
OVERFLOW if, 

avfraof valuf of ovfrflow 



PFR CENT 
OF TOTAL 
11 . 11 
88 . 89 

19.625 



ClJflULATI VF 
PERCENTAGE 
11.11 
ion. on 



HlMIJLATI VF 
R r f 1 A | NP FR 
88.80 
0.00 



SUM OF ARGUMENTS 
314.000 



MULTI PLF 
OF MEAN 
0.057 



DEVIATION 
FROM MFAN 
-1.301 



print the car table statistics. 



TABLE 2 



ENTRIES IN TABLE 
43 



MF AM ARGUMENT 
13.G28 



STANOARP DEVIATION 
13.175 



SUM. OF ARGUMENTS 
586.000 



UPPCR OBSERVED 

L f M T T FRFQUFMCY 

1 3 

OVERFLOW • 40 

AVFRAOF VALUF OF OVERFLOW 



PFR CFNT 
OF TOTAL 
G.98 
93.02 



OUMUL ATI V p 

percentage 

6.98 

100.00 



14.650 



CUMULATIVE 

REMAINDER 

93.02 

0.00 



MULTIPLE 
OF MFAN 
0.073 



DEVIATION 
FROM MEAN 
-0.958 
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what is the current contents of the pump? 

THE CURRENT CONTENTS OF THE PUMP IN THE STATION IS 0. 
what is the s of the pump? 

THE CURRENT CONTENTS OF THE PUMP IN THE STATION IS 0. 

(Note that the above questions to the system are equivalent. 
Questions can be asked by specifying the GPSS standard 
numerical attribute ("s") or by specifying that attribute 
in English ("current contents"). ) 

what is the sm of the pump? 

THE MAXIMUM CONTENTS OF THE PUMP IN THE STATION IS 1. 



what is the sa of the pump? 



THE AVERAGE CONTENTS OF THE PUMP IN THE STATION IS 0. 659. 



what is the st at the pump? 

THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 
10.729. 



what is the qz at the pump? 

THE NUMBER OF ZERO ENTRIES AT THE PUMP IN THE STATION 
IS 38. 



what is the qx at the pump? 

THE AVERAGE WAITING TIME EXCLUDING ZERO ENTRIES AT THE 
PUMP IN THE STATION IS 12. 714. 



what is the maximum line at the pump? 

THE MAXIMUM LINE LENGTH AT THE PUMP IN THE STATION IS 2. 



what is the current line at the pump? 

THE CURRENT LINE LENGTH AT THE PUMP IN THE STATION IS 0. 
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what is the mean transit time of the cars? 

THE MEAN TRANSIT TIME OF THE CARS IS 13.627. 
what is the average line at the pump? 

THE AVERAGE LINE LENGTH AT THE PUMP IN THE STATION IS 0. 278. 
what is the average utilization of the pump? 

THE AVERAGE UTILIZATION OF THE PUMP IN THE STATION IS 0. 659. 
what is the average waiting time at the pump? 

THE AVERAGE WAITING TIME AT THE PUMP IN THE STATION IS 
4. 525. 

what is the average time at the pump? 

THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 
10. 729. 



continue the simulation. 

(This is comparable to the encountering of a second or 
subsequent "start" card in GPSS. The simulation pro- 
gresses from the point at which it previously terminated. ) 

SIMULATION TIME IS 1920(RELATIVE), 1920(ABSOLUTE). 

what is the average utilization of the pump? 

THE AVERAGE UTILIZATION OF THE PUMP IN THE STATION IS 
0. 654. 

what is the average waiting time at the pump? 

THE AVERAGE WAITING TIME AT THE PUMP IN THE STATION IS 

6 . 111 . 

what is the average time at the pump? 

THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 
10. 930. 
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(The three questions above were asked in order to show 
the effect the longer simulation time had on the previous 
results . ) 



continue the simulation 



SIMULATION TIME IS 2880(RELATIVE), 2880(ABSOLUTE). 



what is the average utilization of the pump? 



THE AVERAGE UTILIZATION OF THE PUMP IN THE STATION IS 
0. 701. 

what is the average waiting time at the pump? 

THE AVERAGE WAITING TIME AT THE PUMP IN THE STATION IS 
8. 674. 



what is the average time at the pump? 

THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 
12. 242. 



reset and start the simulation. 

(This is similar to using the control cards "res et/start M 
in GPSS. ) 

SIMULATION TIME IS 960(RELATIVE), . 2840(ABSOLUTE). 

what is the average utilization of the pump? 

THE AVERAGE UTILIZATION OF THE PUMP IN THE STATION IS 
0. 768. 



what is the average waiting time at the pump? 

THE AVERAGE WAITING TIME AT THE PUMP IN THE STATION IS 
9.732. 
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what is the average time at the pump? 



THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 

12 . 982 . 



clear and restart the simulation. 

(This is similar to using the control cards "clear /start" 
in GPSS. ) 

SIMULATION TIME IS 960(RELATIVE), 960(ABSOLUTE). 

what is the average utilization of the pump? 

THE AVERAGE UTILIZATION OF THE PUMP IN THE STATION IS 
0. 732. 



what is the average waiting time at the pump? 

THE AVERAGE WAITING TIME AT THE PUMP IN THE STATION IS 

14. 132. 



what is the average time at the pump? 

THE AVERAGE SERVICE TIME AT THE PUMP IN THE STATION IS 
13. 557. 



are there two pumps in the station? 

(This question results in an inspection of the IPD to 
determine the number of pumps. ) 

NO, THERE IS 1 PUMP IN THE STATION. 

the station should have two pumps. 

(This statement results in an updating of the IPD to 
reflect this new information. ) 
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how often does a vehicle arrive? 



THE TIME BETWEEN ARRIVALS OF THE VEHICLES AT THE 
STATION IS NORMALLY DISTRIBUTED WITH A MEAN OF 8 
MINUTES AND A STANDARD DEVIATION OF 2 MINUTES. 



the mean of the time between arrivals should be 3 minutes, 
and the deviation of the time between arrivals should be 1 



minute. 



develop the x vector. 

(Since the IPD has been altered, the X-vector must be 
reinitialized prior to simulation to reflect the changes 
made. In this case a GPSS program printout is not 
desired. ) 

perform the simulation. 

(This results in the simulation being performed with 
the X_vector for the modified problem. ) 

SIMULATION TIME IS 960(RELATIVE) , 960(ABSOLUTE) . 

what is the average utilization of the pumps? 

THE AVERAGE UTILIZATION OF THE PUMPS IN THE STATION IS 

0. 093. 

what is the average waiting time at the pumps? 

THE AVERAGE WAITING TIME AT THE PUMPS IN THE STATION IS 

3.882. 

what is the average time at the pumps? 

THE AVERAGE SERVICE TIME AT THE PUMPS IN THE STATION IS 

10.529. 
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IV. IMPLEMENTATION 



Achieving the objectives of this research required making additions 
and changes to the FORTRAN routines of NLPQ and to the rules and 
declarations processed by those routines. This chapter, which is 
intended to explain the modifications, has been divided into two sections. 
The first section describes the FORTRAN modifications, and the sec- 
ond section describes the rule modifications. 

e 

A. FORTRAN ROUTINES AND MODIFICATIONS 

The main FORTRAN programming effort involved major altera- 
tions and additions to the simulation routine. A few additional modifica- 
tions to NLP were needed, however, in subroutines PRINT, ENCODE, 
and CRSEG. A discussion on each of these four areas follows. 

1. NLP Modifications to Allow X- Vector Read and Write 

This section was added simply for the convenience of the 
user. It gives the user a method for saving the contents of the 
X-vector as a binary file for use in a later terminal session. By doing 
this the user does not have to duplicate his efforts from a previous 
session to arrive at the same point in a given simulation. He need 
only initialize the current X-vector with the contents of the previous 
X-vector. 

The FORTRAN routine for performing this function is shown 
in Appendix A. This routine was added as entry point XRDWR 
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("X-vector Read/Write ") in subroutine PRINT. The routine may be 

invoked at any time as a command to NLPQ. The format for a call to 

XRDWR from the terminal is 

X A READ A f: or X^WRITE^f: , 

where M f M is an integer number of any available file under CP/CMS 

and denotes at least one blank space. As an example, the 

command M X WRITE 3: n would cause the current contents of the 
/ 

X-vector to be written into file FT03F001 as a binary file. 

2. Establishing a Linkage for Communication 

The rule language of NLP utilizes declared ROUTINES to 
communicate with the various FORTRAN subroutines. The appearance 
of a routine name in a rule causes execution of the code for that par- 
ticular routine in subroutine CRSEG ("Create Segment"). To establish 
communication between the rules and the simulation subroutine, there- 
fore, it was necessary to add an additional routine in subroutine 
CRSEG to communicate with each entry point in the simulation sub- 
routine. The FORTRAN code for establishing this communication is 
given in Appendix A. 

The four entries into the simulation subroutine (SIMULT, 
SIMOUT, SETIND, and SPSTAT) will be discussed in detail later in 
this section. At this point it is sufficient to say that SIMULT 
("Simulation") and SIMOUT ("Simulation Output") require no informa- 
tion from the rule segment being processed. In addition, they are 
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called merely to perform their functions and not to return a value. 
Hence, execution of these routines in CRSEG simply results in a 
call to the appropriate entry point in the simulation subroutine. 

Both SETIND ("Set Indicators") and SPSTAT ("Specific 
Statistics"), however, require the value of certain attributes of the 
segment being processed to be passed as arguments. In addition, 
SPSTAT returns a value which must be made available to the segment 
being processed. The additional coding in CRSEG for these two 
routines sets up this communication. 

3. Modified Output Routine 

The output routine in subroutine ENCODE was originally 
implemented to handle only integer half-word values (or integer half- 
word values expressed in parts per thousand) and, therefore, could 
output only integers in the range -32768 to 32767 or real values in 
the range -32.768 to 32.767. This was acceptable when the only output 
from the system was a GPSS program. With the incorporation of the 
simulation routine and the ability to actually run the simulation at 
the terminal and question the system for such items as the mean 
transit time or the average utilization, however, this limitation 
became unacceptable. As a result, the output routine was rewritten 
to handle the magnitude of any integer or real value which could be 
stored in a full-word on the IBM 360. The modified section of 
ENCODE pertaining to the output of numerical values is shown in 
Appendix A. 
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To output a numerical value from an OUTPUT segment, the 
segment must have an attribute (ATTR or @) 14, 15, or 16. An 
integer or real value will then be output in accordance with one of the 
four cases described below. 

Case 1. If an @14 cell is present and the TYPE of the 
cell is M 0 n , then the value in the address (ADDR) field 
is taken to be a half-word integer. 

Case 2. If an @14 cell is present and the TYPE of the 
cell is n l n t then the value of the ADDR field is taken to 
be a pointer to another cell whose value is a full-word 
found in the ADDR and LINK fields. This value is taken 
to be an integer unless a "l" is present in the ATTR 
field of the same cell, in which case the number is 
considered to be floating point. 

Case 3. If an @15 cell is present, then the value in 
the ADDR field is taken to be a real number expressed 
as a half-word integer in parts per thousand. 

Case 4. If an@16 cell is present, then the value in 
the ADDR field is taken to be a pointer to another cell 
whose value is a full-word found in the ADDR and 
LINK fields. This number is evaluated as in Case 2 
above. 

These four methods of handling numerical values in an 
OUTPUT segment are illustrated in Figure 2 # 
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4. The Simulation Subroutine 



As previously mentioned, modifications to the simulation 
subroutine (called SIMULT) constituted the major FORTRAN program- 
ming effort involved in the preparation of this thesis. A complete 
listing of the subroutine is given in Appendix C. This listing is 
preceded by an extensive dictionary of variables (Appendix B) to assist 
the reader in understanding the logic of the program. Even though 
SIMULT is only a subroutine of NLPQ it is quite extensive and requires 
a considerable amount of core. The source deck contains approxi- 
mately 1600 cards, and a region of 250K is needed for compilation 
under OS/360. A slightly larger region is required under CP/CMS. 
Approximately one and one-half minutes of CPU time are required 
for compilation using the FORTRAN G-level compiler, and the 
resultant object module occupies approximately 40K bytes in the 
IBM 360. The subroutine contains four entry points. It can be 
accessed by calls to SIMULT, SIMOUT, SETIND, or SPSTAT. Each 
of these sections will be discussed individually at this point, 
a. SIMULT 

This is the main entry into the simulation subroutine. 

A call to SIMULT is a request to perform the simulation based on the 
information found in the X-vector. For this reason the user must 
ensure that the vector has been properly initialized prior to request- 
ing that a simulation be performed. If the user desires to continue 
the problem from a previous session in which the contents of the 
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vector were saved, he may utilize the X READ feature described 
previously to initialize the current X-vector. If this is not the case, 
however, the user must either tell the system to develop an X-vector, 
or he must tell the system to write a GPSS program, in a manner to 
be described later. The former case will result in only the initializa- 
tion of the X-vector based on the information contained in the IPD. 

No output will be sent to the terminal. The latter case will result in 
both the GPSS program being output and the concurrent initialization 
of the X-vector. It is important to note that even though NLP gives 
the user the capability of saving the internal problem description 
and reinitializing the system with that IPD at a later date, this 
procedure does not reinitialize the X-vector. 

The SIMULT section is basically the simulation routine 
written by R. J. Williams [9]. The basic logic flow and the structure 
of the various statistical entity layouts described by Williams were 
retained in the implementation of the simulation capability to NLPQ. 
The reader interested in following the logical flow of the SIMULT 
listing is advised, therefore, to reference Williams* work for the 
physical layout of the various statistical entities (STORAGE, QUEUE, 
TABLE, etc.). 

Major alterations to the basic flow of transactions 
through the system were made upon incorporation of the simulation 
routine into NLPQ. For example, a reordering of the sequences 
involved for merging transactions of equal priority into their 
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appropriate positions on the current and future events chains was 
needed in the verification phase to ensure that transactions would be 
executed in the same sequence as is done in GPSS. Modifications 
were also required in the procedures associated with determining 
when a status change has occurred within the system. These modifica- 
tions further altered the original flow of transactions since these pro- 
cedures determine at what points in the simulation a scan of the 
current events chain should be continued from its previous position 
and at what points it should be restarted from the beginning of the 
chain. 

Additional modifications included the addition of the 
entire section pertaining to the updating of the system performed 
during the final time interval prior to terminating the simulation. 

This area had been completely omitted in Williams 1 work but was 
needed to ensure agreement between the statistical outputs of GPSS 
and subroutine SIMULT. This section primarily performs an update 
of the STORAGE and QUEUE statistics to reflect the effect of the 
simulation time involved between the time the storage or queue last 
changed status and the time the simulation was actually terminated. 

The ability to output a selected portion of the X-vector was also 
included in this section. 

Further modification to the original simulation sub- 
routine involved the addition of code throughout the routine to avoid 
system interrupts, such as divide checks, when working with empty 
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storages, queues, and tables. Similar modifications were also 
required throughout the routine to ensure that null pointers in the 
various directories were recognized as such and not used as valid 
indices when storing or retrieving information pertaining to the sub- 
scripted X-vector. 

The ability to properly handle clear, reset, and continue 
commands by the user required some alterations to properly re- 
initialize the allocated storages, queues, and tables. In addition the 
procedure for processing the TERMINATE block was modified to 
ensure replacement of the timing loop generate block on the future 
events chain. The use of both an absolute and relative clock during 
the simulation was deleted and replaced by a single clock (absolute) 
for use during the simulation. A base clock is set in the RESET area 
to allow computation of the relative time in the two instances in 
which a relative time is required, i. e. , (1) the "Cl 11 Standard Numer- 

r 

ical Attribute and (2) the message signaling completion of the 
simulation. 

The section for argument evaluation wa s altered and 
expanded somewhat to allow requests for specific statistics (entering 
the routine via SPSTAT) to access that area of code for evaluation 
of the statistical information requested. Most of the remaining 
alterations to the original routine were either minor in nature or 
made simply in the interest of cleaner coding. 
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SIMULT is entered upon a user request to perform the 
simulation or a request to clear, reset, or continue the simulation. 
Initialization of the simulation model is then performed, if necessary, 
based on the type of request made. This initialization is analogous 
to that performed by GPSS upon encountering the control cards 
SIMULATE/START, CLEAR, RESET, and START. The algorithms 
which direct the flow of transactions through the system from this 
point are essentially the same algorithms used in GPSS. Since the 
results produced by the program, as well as the information needed 
to perform the simulation, are contained in the X-vector, execution 
of SIMULT results in an X-vector altered to reflect the current 
status of the simulation. The results, therefore, are readily avail- 
able to be accessed in the event the user requests information con- 
cerning the outcome of the simulation. Assuming no error conditions 
are encountered during the simulation, the only output from a SIMULT 
entry will be a message giving the absolute and relative simulation 
clock times. This message signals completion of the simulation, 
b. SIMOUT 

The "Simulation Output" routine is invoked whenever the 
user requests GPSS-like statistical information. The format of the 
information output by SIMOUT is practically identical to that of GPSS. 
Unlike GPSS, however, SIMOUT has the capability of providing the 
user with (1) an entire statistical printout (including storage and queue 
statistics, tables, savevalues, the current events chain,' and the future 
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events chain), (2) a single block of any of the statistics just mentioned 
(the storage statistics alone, for example), (3) a single line of storage 
or queue statistics for a single table of several tables (for example, the 
queue statistics for a specified stationary entity), or (4) the queue and 
storage statistics for any specified stationary entity (such as a pump). 
For example, a user request to "print the GPSS statistics" will result 
in (1) above. Similarly, "print the storage statistics" will result in 
(2) and "print the pump queue statistics" will result in (3). If the user's 
request is "print the pump statistics", the SIMOUT routine will output 
both the storage and queue statistics for the pump. This will be explain- 
ed more fully later. 

The SIMOUT routine has no access to the current segment 
being processed. Yet the routine needs to know exactly which lines are 
to be output and which are not. This information is obtained from a 
vector called STATSW ("Statistic Switches"). This vector is local to 
the SIMULT routine and, hence, can be accessed by both SIMOUT and 
SETIND. It is the function of SETIND (which will be described later in 
this section) to set the proper statistic switches for SIMOUT. A call to 
SIMOUT, therefore, must always be preceded by a call to SETIND. 

These calls, however, result from the processing of the input text and 
are made from subroutine CRSEG. They are completely transparent to 
the user. 

STATSW (shown in figure 3) is a four element, real*8 
vector. The first three elements contain indicators for storages, 
queues, and tables, respectively. The rightmost 55 bits in each 
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STORAGE 

QUEUE 

TABLE 



(1-55) 

(1-55) 

(1-55) 



Figure 3: Statistic Switch Vector. 



43 



element are used to indicate which storage, queue, or table should be 
output. If the first element in STATSW has a "1" in bit position two 
and zeros elsewhere, for example, the SIMOUT routine would output 
the storage statistics for the stationary entity which has an identifica- 
tion number (IDNO) of two in the IPD. The storage statistics for the 
remaining stationary entities would not be printed. Indicators for 
queues and tables are treated similarly, with the exception that the 

i 

IDNO for the table refers to a mobile entity in the IPD. 

The fourth element in STATSW contains indicators for 
savevalues and the current and future event chains. Only the right- 
most three bits are used. A "l" in the third bit position of element 
four, for example, would indicate the future events chain is to be 
output. The bit pattern shown in figure 3 indicates that storage and 
queue statistics for the stationary entity having an IDNO equal to 3 is 
to be output. No other statistics are to be printed. 

Upon entry into SIMOUT a call is immediately made to 
subroutine GBITS ("Get Bits"). This routine returns the value of bits 
1 through 55 of the first STATSW element. A zero value indicates no 
storage statistics are to be output. In this case a branch is made 
around the "output storage" area to the "output queue" area. If the 
value returned is not zero, however, then one or more lines of storage 
statistics are to be output. SIMOUT, therefore, outputs the storage 
headings and then begins a search to determine which of the storages 
are to be output. This is done by successive calls to GBITS to obtain 
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the value of bit positions 1, 2, 3, ... (up to the number of storages 
being maintained in the system). If the individual bit value is "0", 
that storage is bypassed. If the bit value is "1", however, SIMOUT 
utilizes the current information in the X-vector for that particular 
storage to calculate and output the statistics required. When this has 
been done for each storage, the program falls through to the output 
queue" area. 

The procedure for determining which queues and tables 
are to be output, if any, is the same as that for storages. Upon 
falling through to the "output savevalues" area, however, only one 
call to GBITS is made to determine if bit position one of STATSW(4) 
is on. If this is the case, all of the savevalues are listed with their 
respective values. 

The procedure for determining whether the current and 
future events chains (bit positions two and three) are to be output is 
essentially the same as that for savevalues. The noticeable difference 
in output of the chains is that both are handled in the same area in 

order to avoid duplication of code. 

Upon completion of the statistical outputs, SIMOUT zeros 

all four elements of STATSW prior to returning to CRSEG. This is 
required to prevent unwanted statistics from being printed if the user 
later requests further information requiring another call to SIMOUT. 
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c. SETIND 



As mentioned previously, it is the function of SETIND 
( n Set Indicators’ 1 ) to set the desired bits in the STATSW vector to 
enable SIMOUT to output the proper statistics. It is a function of the 
decoding rules to determine the content of the English request and 
issue a call to SETIND, telling the routine which bits are to be set. 

Like all calls to the simulation routine, this call is made via sub- 
routine CRSEG. 

SETIND must receive two parameters from the calling 
segment. These parameters are (1) the row in STATSW which is 
affected and (2) the bit position (or positions) within that row which is 
to be set. Since a single request for statistics by the user may involve 
the setting of a single row or multiple rows and may involve the set- 
ting of a single bit or all the bits within a given row, this capability 
has been included in SETIND in order that these multiple settings 
might be performed by a single call to SETIND. This is accomplished . 
by the manner in which SETIND handles the calling arguments. If the 
requested row is in the range 1 through 4, the routine assumes the 
row specified is to be set. If the calling argument for the row is 5, 
however, the routine assumes that the bit (or bits) specified are to be 
set in both rows 1 and 2 of STATSW. This condition is common to a 
request for statistics of stationary entities (which have both storage 
and queue statistics). A calling argument of 6 specifies that all the 
bits of each element are to be set. This corresponds to a total 
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GPSS-like printout. If the second calling argument (the bit position 
to be set) is in the range 1 through 55, that single bit position is set. 
If the second argument is greater than 55, however, then all of the 
bits of the specified row are turned on. This condition is common 
when issuing a request for statistics without specifying an associated 
stationary or mobile entity. 

The actual setting of the bit positions is not necessarily 
performed by SETIND. SETIND merely analyzes the request to 
determine which bits are to be set. If an entire row in STATSW is 
to be turned on, SETIND will perform the function. If only a single 
bit is to be turned on, however, SETIND issues a call to subroutine 
PBITS ("Put Bits") to set the desired bit position. Once the desired 
bits have been set, SETIND returns to the calling routine, CRSEG. 
d. SPSTAT 

A user request for a single item of statistical informa- 
tion is processed by the rules in a manner similar to any other ques- 
tion to the system. The value of most items of statistical interest, 
however, is not available in the IPD. When it is determined in the 
processing of the text that one of these values is being requested, 
attributes are set in the segment being processed to designate the 
type of value being requested and the IDNO of the entity for which 
the value is to be computed. A call is then made through CRSEG 
to SPSTAT ("Specific Statistics") passing these attributes as 
parameters . 
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SPSTAT sets an initial default value of zero to be returned 



in the event an error condition occurs (such as a request for statistical 
information prior to performing the simulation). The routine then sets 
an entry point flag and branches into the SIMULT argument evaluation 
section to compute the desired statistic. As previously mentioned, 
this section in SIMULT was modified to be able to process inputs from 
both entry points. A SIMULT entry into this section causes all real 
argument values to be truncated to integer values. An SPSTAT entry, 
however, requires that real values be retained and returned as float- 
ing point. 

Upon completion of argument evaluation, the entry point 
flag directs the logical flow back to SPSTAT. The bit pattern of the 
value is set into an integer word and a flag is set to specify whether 
the value being returned should be interpreted as an integer or a 
decimal result. The requested value and the flag are then passed 
back to CRSEG which inserts the information into the current segment 
record to be output later in the encoded text. 

B. RULE ADDITIONS AND MODIFICATIONS 

Integration of the simulation routine and the ability to query the 
system as to the results of the simulation required several modifica- 
tions and additions to the existing rule modules. Expansion of the 
concept structure by the addition of named record definitions was 
also necessary to allow questions relating to the simulation results. 
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The changes and additions made are shown in Appendix D and are 
discussed in the following sections, 

1 . Named Records 

Several named records in the form of English words with 
their associated part-of-speech (PS) attributes were added to facilitate 
recognition of these words by the system. The GPSS entities, storage, 
queue, and table, were also declared and assigned numerical codes 
which correspond to the position of their respective indicators in the 
STATSW vector previously described. A superset relation is also 
established to identify these words as elements of the set 'GPSSENTY' 
("GPSS entity"). 

The remaining named record definitions identify the various 
GPSS "Standard Numerical Attributes" (SNA's) as members of the 
set 'GPSSATTR' ("GPSS attribute"). Each of these records also 
contains an SNACODE attribute with an associated value. This value 
is passed to the SPSTAT routine by the rules in those instances where 
a specific statistic is desired. The CHARS attribute contains the 
SNA name in character form to be used in encoding responses to the 
user's questions. 

2. Simulation Control Commands 

To permit interactive control of the simulation with regard 
to the various modes of operation, several semological decoding 
rules were required. These rules are basically "key word" rules 
which serve to test the input string for the presence of one or more 
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key words. Rules of this type have a left segment type of KWDSENT. 
The condition specifications of these rules indicate the key words 
necessary for rule execution. For example, the presence of the key 
word "perform”, "simulate", or "run", in the user ! s request results 
in execution of the simulation routine in the SIMULATE/START mode. 
Thus the commands; "Perform the simulation. ", "Simulate the 
system. ", or just "run. ", are equivalent and result in the simulation 

i 

being run. 

The keywords "reset", "clear", and "continue" are handled 
in a similar manner. Thus by using commands such as "Reset and 
start the simulation. ", "Clear the model. ", or " Continue the 
simulation. ", the user can control the mode of the simulation. The 
key word rules which handle these cases set the mode indicator of 
the X-vector (X(l)) to the appropriate value and reset the termination 
count. In the present application, the termination count is set to 1 
since the GPSS/X- VECTOR rules use a "timing" transaction to 
terminate the simulation. The simulation is automatically restarted 
once these modifications have been made. 

Several other functions are also performed by the key word 
rules. The key words "gpss" or "vector" occurring in the input 
command result in the initialization of the X-vector from the IPD. 

If only the key word "vector" is present, for example "Develop the 
x vector . ", the GPSS program will be suppressed. If "gpss" is 
present, both the GPSS program and the X-vector initialization will 
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result. The presence of the key words "print" and "gpss" combined 
with the absence of the key word "program" produce a complete GPSS 
statistical listing. Thus, "Print the gpss statistics. " and similar 
constructions produce output with the complete results of the 
simulation. The key words "print", "current", and "events " appear- 
ing in the input text result in a printout of the transactions on the 
current events chain. Omission of, or substitution for, the key word 
"current" produces a listing of the future events chain. Combination 
of key words which satisfy more than one rule (e. g. , "Reset and run. ") 
will result in execution of the first applicable rule based on their 
physical order as shown in Appendix D. 

3. Producing Selective Simulation Results 



opportunity to request certain portions of the GPSS-like statistical 
printout. The general format for commands of this type is: 



Thus commands of the form "Print storage statistics. " result in 
that portion of the statistical printout related to the GPSS entity 
specified; in this case, the storages. The GPSS entities which can 
presently be employed are "storage", "queue", and "table". The 
command "Print pump statistics. " satisfies the general format and 
produces all statistical output related to the stationary entity "pump". 

In the present application, the statistical output produced for stationary 



Several rules were added to NLPQ to give the user the 



PRINT 




STATISTICS. 
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entities is the appropriate line of storage and queue statistics, with 
their respective headings. Substitution of a mobile entity in the same 
form (e. g. , "Print truck statistics.") produces table statistics for 
that mobile entity. 

Further selectivity can be obtained by supplying more optional 
information as in the command "Print the pump queue statistics. " In 
this- instance only the line of output associated with the queue at the 
pump will be printed. The corresponding command for mobile entities 
(e. g. , "Print the truck table statistics. ") is equivalent to the earlier 
mobile entity command and also results in a table. Care must be 

taken when utilizing this form to ensure that the selection of entities 

\ 

is compatible. Stationary entities must be used in context with the 
GPSS entity "storage" or "queue". Mobile entities require the use 
of the GPSS "table" entity. Incorrect sequences may result in alter- 
nate statistics or none at all. 

The rules for processing these "print commands" are includ- 
ed in the portion of the "Lexology for Decoding English" shown in 
Appendix D. The processing is based on the appearance of noun 
phrases which are elements of either the set 'GPSSENTV or 
'ENTITY'. A noun phrase (NOUNPH) which is in the set (denoted by $) 
'GPSSENTY', such as "the storage", results in a STATPH segment 
containing the appropriate code for that entity in attribute eight. 
Occurrence of the STATPH segment in the context "print STATPH 
statistics. " results in the creation of a PRINTPH segment which 
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produces the appropriate calls to SETIND and SIMOUT to output the 
block or line of statistics. The calling parameters used for SETIND 
are the values of attributes eight and nine of the PRINTPH segment. 
These attributes and values are copied directly from the STATPH 
segment. 

STATPH segments are also produced by noun phrases in the 
sets 'STATENTY 1 ("stationary entity") and ’MOBENTY 1 ("mobile 
entity"). These instances, such as "the pump" or "the car" in the 
proper context result in the production of the single lines of output 
corresponding to stationary entities or the table statistics for mobile 
entities. A series of two noun phrases, the first of which is in the 
set ■ENTITY 1 (either stationary or mobile) followed by a noun phrase 
in the set 'GPSSENTY' (storage, queue, or table) also results in the 
creation of a STATPH. In this case, attribute eight is set to the 
code of the GPSS entity (accessed via the second noun phrase), and 
attribute nine is set to the identification number of the entity (via the 
first noun phrase). These parameters are then used in the SETIND 
routine to indicate the pertinent line of statistics to be produced by 
SIMOUT. A comparison of the rules and the general format described 
earlier provides an insight into the way in which these commands are 
presently handled. 

4. Interrogating the Simulation Results 

Additional rules were also incorporated to allow the user to 
ask questions about specific results of the simulation. With this added 
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capability, total, partial, or individual statistics are readily avail- 
able to the user. Presently acceptable questions are of the form: 

C ^ f SNA name! I OF*) j ) [Stationary Entity) 

WHAT IS JtheJ |g NA J | at J |THE| [_^ obile Entity J 

The SNA's currently in use are those which correspond with the 
individual statistical elements produced in the GPSS-like printout. 
They are listed in the named record definitions in Appendix D 
(beginning with 'SC'). The character string attribute (CHARS) of 
each of these records serves as a natural language "SNA name 
which can also be used in most instances. In utilizing the question 
form above, the user again must ensure compatability between the 
SNA (or SNA name) and the type of entity statistic desired. 

The two questions, "What is the_SR of the pump ? " and 
"What is the average utilization of the pump ? ", are equivalent and 
produce a response with the appropriate number. The same is true 
of the questions "What is the TB of the trucks? " and "What is the 
mean transit time of the trucks? ". The only SNA's available 
presently for the mobile entities are TB, TC, and TD, which are 
"mean transit time", "number of entries", and "standard deviation". 
Storage and queue individual results may be obtained by using the 
SNA's SC, SM, SR, SA, S, R, ST, Q, QA, QM, QC, QZ, QT, QX, 
or QP with their associated stationary entity. The English name of 
each of these SNA's is contained in the CHARS attribute of the corres 
ponding named record in Appendix D. 
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The rules for recognizing SNA names are also contained in 
the decoding lexology. Individual rules exist for each allowable SNA 
name. For example, in the question, "What is the current line at 
the pump? ", "current line" results in a noun phrase segment in the 
set f Q r . This noun phrase segment is later processed by an encoding 
rule (QUEST2) which tests for this set relation. Satisfaction of the 
rule conditions result in a sentence segment with attribute eight con- 
taining the SNACODE of the desired statistic and attribute nine contain- 
ing the IDNO of the mobile or stationary entity. Using the values of 
attributes eight and nine as calling parameters, SPSTAT is called to 
return the value of the desired statistic. The value returned is then 
used in the response to the user. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The thesis objective has been met. The simulation routine and 
associated rules have been integrated into the existing NLPQ system 
to produce an initial version of an interactive simulation system. 

With this feature, the NLPQ user can perform, and control the mode 
of, the simulation for his specific problem by using natural language 
commands. The results of the simulation may be requested in 
several ways. A complete GPSS-like printout is available or the user 
may select those portions of the printout which are of interest in his 
specific problem. Blocks of statistics (storage, queue, table, etc. ) 
or individual lines of the statistical output are readily accessible by 
the user in the latter instance. Specific statistical results may also 
be requested in a question-answer fashion. Using these additional 
features, the user can solve queuing system problems in an inter- 
active manner through natural language dialogue with the system. 

Recommendations for further research include expansion of the 
present question handling abilities to allow further interrogation of 
the simulation results in a less stringent manner. The ability to 
detect incompatibilities in input questions and commands and produce 
meaningful error analyses would enhance the present interaction with 
the user. Extension of the present features of the simulation routine 



56 



to include a larger subset of those functions performed by GPSS, 
combined with the necessary rule language additions, would provide 
the user- with greater power in solving more complex problems. 
Further enlargement of the present rule modules to permit inter- 
active ability in specifying table limits, run times, and other GPSS 
attributes, would enable greater specification and control of the 
simulation by the user. A means of handling situations in which 
aggregate statistics are desired (for example, the station statistics 
should reflect the aggregate statistics for the pump or pumps in the 
station), would also be of value in improving the usefulness of the 
system. 
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APPENDIX A 



*** NLP MODIFICATION TO ALLOW *** 
X VECTOR READ OR WR ITE 



C XRDWR 

ENTRY XRDWR 
501 J=J+1 

IF (COL (J) .EQ. BLANK) GO TO 501 

IF CCCLI J) .EQ. COLON) CALL ERRORA ( J , 9 , £5 50 ) 

CALL COLECT (NAME) 

511 J=J+1 

IF ICOLIJ) . EQ. BLANK ) GO TO 511 

IF (COLI J) .EQ. COLON) CALL ERRORA t J ,9 ,£550 ) 

CALL CCNVRT ( NUM ) 

IF (NUM. LT.l.OR.NUM.GT .14) CALL ERRORA ( J , 10, £550 ) 
NUM4=NUM 

IF (NAME.EO.XRD) GO TO 531 
IF (NAME. EO.XWR) GO TO 541 
CALL ERRORAI J, 9, £550) 

531 REWIND NUM4 

READ ( NUM4 ) { X ( I ) , I =1 , MAXX ) 

RETURN 

541 REWIND NUM4 

WRITE { NUM4 ) { X ( I ) , I =1 , MAXX ) 

550 RETURN 
END 



*** ADDITICNAL ROUTINES *** 



C 

3210 

C 

3220 

C 

3230 



C 

3240 



CALL S I MULT { TR6 , CUT6 , RT ERM , WT ERM ) 
GO TO 1200 

CALL SIMCUTI0UT6) 

GO TO 1200 

R0W=HVAL(A8,SEGMNT) 

BIT=HVAL(A9,SEGMNT) 

CALL S ET I ND ( ROW ♦ E I T ) 

GO TO 1200 



SIMULT 



S IMOUT 



SETIND 



SPSTAT 



SNA=HVAL { A8 , SEGMNT ) 

IDT=HVAL(A9 .SEGMNT) 

CALL SPSTATISNA, IDTtVAL, IORD) 

LX=L0C(A9, SEGMNT ) 

CEL=CELL (LX) 

TYPE=1 

ADDR = NEWCEL ( CONST R( ZERO, I ORD , PV AL ( 1 ) , PV AL ( 2 ) ) ) 
CELL { L X ) =C EL 
GO TO 1200 
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*** MODIFIED OUTPUT ROUTINE *** 



C PROCESS 'OUTPUT' 

C 

201 N = HVAL(AlltSEGMNT) 

IF (N.GT.O) CALL OUTCHR ( DZERO) 

N = N - 1 

IF (N.GT.O) CALL SKIP(N) 

KOL = HVAUA12 , SEGMNT) 

IF (KOL.GT.O) CALL SETJJ(KOL) 

NXTWRD = 0 

FWRD = LOC (A13 , SEGMNT ) 

IF (FWRD.EO.O) GO TO 231 
2 21 WORD = DVALUEI FWRD , NXTWRD) 

DO 225 1=8,64, 8 

CHR = DOR(DLS (DRS(W0RD,64-I ),56),DZB) 

225 IF ( CHR.NE. DBLANK) CALL OUTCHR(CHR) 

I f= (NXTWRD. NE.O) GO TO 221 
231 LOCC = LOC (A14 , SEGMNT ) 

IF (LOCC.EO.O) GO TO 241 
DCM=0 

CEL=CELL( LOCC) 

IF (TYPE.EO.l) GO TO 253 

RN=ADDR 

GO TO 257 

241 LOCC = LOC (A15, SEGMNT) 

IF (LOCC.EO.O) GO TO 251 
DCM= 1 

CEL=CELL (LOCC ) 

RN=ADDR/ 1000 . 

GO TO 257 

251 LOCC = LOC( A16, SEGMNT ) 

IF (LOCC.EO.O) GO TO 131 
CEL=CELL (LOCC) 

253 CEL=CELL { ADDR ) 

DCM=ATTR 

IF (DCM.EO.l) GO TO 255 
RN= 14 
GO TO 257 
255 RN=R4 

257 IF (RN.GE.O.) GO TO 261 
RN=-RN 

259 CALL OUTCHR(DMINUS) 

261 IF (RN.GT.1.0E10) GO TO 298 

s w=o 

DO 265 1=1, 13 
11=10-1 

IF ( DCM.EO.O. AND. I I .EO.-l ) GO TO 131 
M=RN/10 .*=<=1 1+0 .0005 

IF ( SW.EQ. D.AND.M.EQ.O.AND.I I.GT.O) GO TO 265 
SW= 1 

IF ( DCM.EO.l .AND. I I .EO .-1) CALL OUTCHR ( DECPNT ) 

CALL OUTCHR(DIGIT( M+l) ) 

RN=RN-M*10.**I I 
265 CONTINUE 
GO TO 131 

298 WRITE ( OUT 6 , 29 9 ) RN 

299 FORMAT (' NUMBER EXCEEDS ''OUTPUT" LIMIT. VALUE IS ' 
1 , E 13. 7) 

CALL OUT CHR ( DS T AR ) 

CALL OUTCHR (DSTAR) 

CALL OUTCHR (DSTAR) 

GO TO 131 
END 
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APPENDIX B 





DICTIONARY OF VARIABLES USED 
IN THE SIMULATION ROUTINE 


ALTBLO 


Alternate block number for TEST or GATE 
routines. 


AVAIL 


Amount of storage available in a given storage 
entity. 


BASE 


Value to which spread will be added in 
GENERATE and ADVANCE blocks to determine 
departure time from the FEC. 


BITTS 


Value of requested bit pattern returned from 
call to GBITS. 


BLOCK 


Pointer to word preceding block directory in 
X-vector (BLOCK=X(19)). 


BYTE (*) 


Logical^ 1 temporary variable 
(BYTE( 1 ) = DW ORD( 1 )). 


CBLO 


Pointer to current block being processed. 


CKPNT 


Check point used for following traces when 
debugging. 


CLOCK 


Absolute clock time (CLOCK=X( 1 1)). 


CLOCKB 


Base clock time needed to compute relative 
clock time after RESET. 


CLOCKR 


Relative clock time. 


COMVAL 


Comparison value used in SELECT blocks. 


CQUE 


Pointer to current queue being processed. 


CREATE 


Time differential between current clock and 
time transaction is due to leave the FEC, 
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CSTO 


Pointer to current storage being processed. 


CTRA 


Pointer to current transaction being processed. 


DELAY 


Time delay caused by ADVANCE block. 


DONES 


Real *8 mask consisting of all ones. 


DTIME 


Real*8 time differential between clock and time 
queue or storage last changed status. 


DW* 


See DWORD(*). 


DWORD)*) 


Real*8 temporary variable (DW* = DWORD(*)). 


EPT 


Variable which flags entry point at which 
SIMULT was entered. 


ERR 


Flag set when an error is encountered 
(ERR=X(30)). 


ERRORZ 


Subroutine called to output error message and 
set ERR flag. 


FB 


Variable which specifies which bits to set when 
calling SETIND or the first bit to set when calling 
PBITS . 


FFW* 


See FFWORD(*). 


FFWORD)*) 


Real*4 temporary variable 

(FFW*=FFWORD(*); FFWORDf l)=DWORD( 1 )) . 


FLAG 


Temporary variable used as a logical flag 
when needed. 


FOS(*) 


Stack for evaluating floating point arguments 
(FOS(*)=OS(*)). 


FREWDT 


Width of frequency interval in a table entity. 


FRNG 


Number of frequency intervals in a table entity. 


FTCEC 


Pointer to the first transaction on the Current 
Events Chain ( FTCEC=X(26)) . 
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FTFEC 


Pointer to the first transaction on the Future 
Events Chain ( FTFEC=X(2 8) ). 


FTRA 


Pointer to the first transaction on the list of 
unused transactions. 


FUNCT 


Pointer to the word preceding the function 
directory in the X-vector (FUNCT=X(17)). 


FW * 


See FWORD(*). 


FWORD(*) 


Integer'*'4 temporary variable 
(FW*=FWORD(*); FWORD( l) = DWORD( 1)). 


i 

GBITS 


Function which returns value of bits specified. 


GIVEUP 


Number of storage units being made available. 


HIGH- 


Ending entity number to be used in SELECT 
block search. 


HW* 


See HWORD(*). 


HWORDf*) 


Integer *2 temporary variable 
(HW*=HWORD(*); HWORD(l)=DWORD( 1)). 


IDT 


Identification number of entity to be used when 
calling SPSTAT. 


IMAX 


Looping limit for outputing the OS stack when 
tracing argument evaluations. 


INST 


Pointer to current block element being 
processed. 


INTDEC(*) 


Indicates whether corresponding OS/FOS value 
is integer or floating point. 


IORD 


Indicates whether value returned by SPSTAT 
is integer or decimal. 


JBL 


Number of the block whose routine is being 
executed. 


KDEX 


Index used to keep count of number of arguments 
specified for a block. 



62 



KMODE 


Special mode variable used to determine 
how a function is to be evaluated. 


KOUNT 


Counter used as index into function entity 
tables. 


KXVAL 


Integer*4 value of independent variable 
used in function computation. 


KYVAL 


Integer*4 value of dependent variable used 
in function computation. 


LB 


Specifies the last bit to be set in a call to 
PBITS. 


LOW 


Beginning entity number to be used in 
SELECT block search. 


LTCEC 


Pointer to the last transaction on the 
Current Events Chain (LTCEC=X(27)) . 


LTFEC 


Pointer to the last transaction on the 
Future Events Chain (LTFEC=X(29)). 


MAXCTS 


Maximum contents of a storage or queue. 


MAXX 


Maximum number of X-vector elements. 


MDFR 


Modifier used in TEST, SELECT, GATE, and 
ASSIGN blocks. 


MODE 


Indicates whether simulation will begin in 
Simulate /Start, Reset, Clear, or 
Start mode (MODE=X(l)). 


MRNG 


Number of points in a function. 


NBLO 


Number of blocks (NBLO=X( 18)). 


NEWBDT 


Block departure time of transaction being 
merged into the FEC. 


NEWPRI 


Priority of transaction being merged into the 
CEC. 
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NEXBDT 


Block departure time of first transaction 
on the FEC. 


NEXBLO 


Next block number transaction is to enter. 


NFUNCT 


Number of functions (NFUNCT=X( 16)). 


NOUPD 


Indicator used to merge a transaction into the 
CEC without updating the clock. 


NPAR 


Number of transaction parameters 
(NPAR=X(24)). 


NQUE 


Number of queues (NQUE=X( 14)) . 


NSAV 


Number of savevalues (NSAV=X(22)). 


NSTO 


Number of storages (NSTO=X( 12 )). 


NTAB 


Number of tables (NTAB=X(20)). 


NVAR 


Number of GPSS-type variables 
(NVAR=X(3 1 )). 


OLDBDT 


Block departure time of transaction being 
checked in the FEC. 


OLDPRI 


Priority of transaction being checked in the 
CEC. 


OPCODE 


Operation code indicating type of block. 


OS(*) 


.Stack for evaluating integer arguments 

(OS(*) = FOS(*)). 


OUTX 


Specifies optional output file for X-vector 
(defaults to 0 (no output)). 


OUT6 


Specifies optional output file for traces and 
SIMOUT output (defaults to 6 (terminal)). 


PBITS 


Subroutine called to turn on specified bits in 
an indicator. 


PCONTS 


Present contents of a queue. 
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PNTA 


Pointer to last entry on return address' 
stack. 


PNTB 


Pointer to last evaluated argument on 
OS/FOS stack. 


PNTF 


Pointer to the first element of the frequency 
intervals for a table entity. 


PNTS 


Pointer to the first element of the slope values 
of the current function entity. 


PNTT 

i 


Pointer to the transaction being compared 
with the current transaction. 


PNTX 


Pointer to the first element of the independent 
variable values of the current function entity. 


PNTY 


Pointer to the first element of the function 
values of the current function entity. 


PR(*) 


Temporary vector to save integer values in 
the OS stack for output while tracing. 


QUE 


Pointer to the word preceding the queue 
directory in the X-vector (QUE=X(15)). 


RAS 


Return address stack used to temporarily 
hold the location of the next arguments to 
be evaluated. 


RN 


Random number. 


ROW 


Argument to SETIND which specifies which 
status switch row is being set. 


RR* 


Temporary real'-M variables used primarily 
to save values for immediate output. 


RSLOPE 


Slope used in function evaluation. 


RTERM 


Specifies optional input file for entering 
optional data (defaults to 5 (terminal)). 


RVAL 


Value to be returned by SPSTAT if real*4 
result is required. 
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RX(*) 


Variable used to manipulate real' : '4 values in 
the X-vector (RX(*)=X(*)). 


RXVAL 


Real*4 value of the independent variable used 
in function computation. 


RYVAL 


Real*4 value of the dependent variable used 
in function computation. 


SAVE 


Pointer to the word preceding savevalue 
locations in the X-vector (SAVE=X(23 )). 


SC FLAG 


Status change flag. Restarts scan at the 
beginning of the CEC when on. 


SE 


Pointer to the top of the pushdown chain for 
storage empty condition. 


SEED(*) 


Seed used in random number generation 
(SEED(1-8)=X(3-10)). 


SETIND 


Entry point. Called to set status switch 
indicators . 


SF 


Pointer to top of pushdown chain for storage 
full condition. 


SIMOUT 


Entry point. Called to output GPSS-like 
statistics . 


SNA 


Specifies the Standard Numerical Attribute 
to be evaluated on call to SPSTAT. 


SNE 


Pointer to the top of the pushdown chain for 
storage not empty condition. 


SNF 


Pointer to the top of the pushdown chain for 
storage not full condition. 


SPSTAT 


Entry point. Called to output special 
statistics (individual SNA’s). 


SRED 


Pointer to the top of the pushdown chain 
for storage reduction in contents condition. 
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SSIND 

STATSW(*) 

STO 

TAB 

TASGN 

TBLNO 

TBLO 

% 

TCNT 

TCTRA 

TEMP* 

TEST 

TPAR 

TRACE 

TRAN 

TRARG 

TRBLO 

TRCNT 



Scan status indicator. True if the transaction 
is in an active status on the CEC. 

Vector containing status switches which 
indicate which statistics are to be output 
by the SIMOUT routine. 

Pointer to the word preceding the storage 
directory in the X-vector (STO=X(13)). 

Pointer to the word preceding the table 
directory in the X-vector (TAB = X(2l)). 

The value being assigned in an ASSIGN block. 

Number of table being tabulated. 

Temporary variable used to save the pointer 
to the current block. 

Temporary location for saved termination 
count in TERMINATE block. 

Temporary variable for saving the pointer 
to the current transaction. 

Temporary integer*4 variable. 

Value of pointers associated with the delay 
chains . 

Parameter being assigned a value in an 
. ASSIGN block. 

Logical switch for tracing simulation mode. 

Pointer to the word preceding the transaction 
entities (TRAN=X(25 )) . 

Logical switch for tracing argument 
eva luation. 

Logical switch for tracing block routines. 

Indicates the transaction number being assigned 
a new transaction. 
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TRMCNT 


Termination count (TRMCNT=X(2)), 


TRSCN 


Logical switch for tracing scan procedure. 


TRSIZE 


Indicates the size of the X-vector printout. 


TRS1 


Logical switch for tracing chain manipulations. 


TRTN 


Temporary variable used to save the ZRTN 
value . 


TRUPD 


Logical switch for tracing update clock 
procedure. 


TR6 


Logical trace enable switch. 


TTRA 


Temporary variable used to save the pointer 
to the current transaction. 


TYPARG 


Indicates whether the search argument of a 
function entity is integer or floating point. 


TYPX 


Indicates whether the independent variable 
of a function entity is integer or floating 
point. 


TYPY 


Indicates whether the function values of a 
function entity are integer or floating point. 


ULLI 


Indicates the upper limit of the lowest frequency 
interval for a table entity. 


UNITS 


Number of units entering or leaving a queue. 


USED 


Current contents of a storage entity. 


VAL 


Value to be returned by SPSTAT as an 
integer*4. 


VAR 


Pointer to the word preceding the variable 
directory in the X-vector (VAR = X(32)). 


WANT 


Number of storage units required. 


WRTN 


Temporary variable used to save the ZRTN 
value. 



68 



WTERM 


Specifies optional output file 
(defaults to 6 (terminal)). 


WTFAC 


Weighting factor utilized in TABULATE block. 
Defaults to one if not specified. 


X(*) 


Vector used for holding all storages, queues, 
tables, directories, etc. 


ZRTN 


Holds statement number for use in returning 
from various routines. 
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